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AM ANALYSIS OF LANDSAT MSS SCm*TO-SGBNE REGISmTIOli ACCUitACY 

I. IMTItODUCnOil 

This p«p«r dfttcriliea •fforca hy NASA's Nstioosl 3psc& iMhnology Lsbors- 
toriss (N^TL) i Earth Issourccs Lsborstory (niL) and U9>A's Statistical Rape rt ins 
Sarvlca (^S) to analysa tha accuracy of scana-to-scans ragistration of Landsat 
Nultlspactral Scatmar (NS8) data, twalva araas. raprasttatlng four aajor typas 
of land cover categories » ware salactad for this study. Landsat MSS data ware 
nerged by both ERL and SRS. ERL used the Earth Resources Laboratory Applications 
Software (ELi^) systm to parfora the registration, ^lle SRS used the EDITOR 
software system* 

The specific objectives of this study irere: 

1. To deteraine the registration accuracies of the data sets nerged 
by ERL and SRS. 

2. To detemlne whether the registration techniques used were equally 
well suited for the different land cover categories studied. 

3. To deteraine whether one registration technique produced significantly 
better registration. 

4. To detemlne whether a significant relationship exists between the 
root nean square (ms) errors of control points used ic scene-to- 
scene registration and the accuracies deterelned by this study. 
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II. LAMDSAT SATA 


One of the goals of this stt^y ifas to deternine the suitability of the 
ELAS and EDITOR registrat^ techniques for a variety of land cover categories. 
The general surface cover categories chosen were agriculturei rangeland* forest* 
and mixed. The mixed category comprises areas with both forest and agriculture. 

Three test sites were selected for each of the four land cover categories* 
for a total of 12 test sites. Two Landsat scenes were used by ERL for each site. 
Table 1 lists each Landsat scene end its land use category. 

To minimize tape mailing between SRS and ERL* two sites from each of the 
four land cover categories were registered by SRS and the merged data were 
mailed to ERL for c(xnparison of the two registration techniques. The two re- 
gistration techniques were coaq>ared for these eight sites. 


2 


TABLE 1. LANDSAT SCENES AND UND COVER CATEGORIES 


LAND COVER 
CATEGORIES 

LOCATION AND LANDSAT SCENE 



MISSOURI mmr 

S. DAKOTA (32/29) 

ARIZONA (HO/57) 

AGRICULTURE 

21654-16100 

30435-16165 

21640-16313 

21676-16321 

21540-17174 

21630-11203 


KANSAS imm 

IDAHO ifs/m 

OKLAHOMA (30/36) 

RANGELAND 

21980-16270 

30653-16254 

22029-17405 

30558-17404 

22016-16283 

21566-16204 


COLORADO (38/34) 

$, CARQLlNA-il8/2fil 

IDAHO (46Z2Z1 

FOREST 

30517-17141 

21718-17095 

30515-15201 

30623-15181 

21654-17511 

30579-17561 


KENTUCKY (23/341 

N. CAROLINA (17/36) 

MISSISSIPPI (23/39) 

MIXED (FOREST 
AND AGRICULTURE) 

21559-15392 

30556-15474 

22057-15142 

21967-15125 

a667-15442 

n739-15A61 




(XX/XX) * PATH/ROW NUMBERS. 



III. ACCURACY ASSESSMENT HBmimjOGt 

For Mch of th« test altos analysad, tha Landaat scanas iiara aargad by TO. 
using tha ELAS soft%»ara. SRS usad tha S)ITOR softwsra to Mrga sight of the 
sites. Tha accuracy of the merged data sets was then assess^. 

The accuracy assessment procedure consisted of generating gray-scale 
electrostatic plots of the data and selecting teat points using a digitiser. In 
this procedure two gray-scale plots were made of two bands of the base Landaat 
scene for approximately 16 different subscene areas. A third plot wa^' made of 
one band of the merged Landaat scene. Then, for each subscene area, all three 
plots were mounted on a digitizer and matching test points were selected fr<mi 
each plot. Ta the cases of Missouri, Kansas, Colorado, and Kentucky, a fourth 
plot was made for each subscene area of one band of the data merged by SRS. 

Then all four plots were mounted on the digitizer at one time. Thus, these 
four comparisons were made on identical points. The foUr remaining cos|>arl8on8 
are based on independently selected test points. 

The points selected from the two plots of two bands of the base Landsat 
scene were used to estimate the random effect of a husan selecting the same test 
points in this manner. The points selected from the plots of the merged Landsat 
data, when compared to the test points of the base Landsat scene, measured the 
random human error in addition to the registration error. A statistical proce- 
dure was used to separate the two sources of error and estimate the portion 
attributable to misregistration by the registration software. 
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IV. STATISTICAL WHm. 

In the aodeJ , . repreeente the set of x end y coordi* 

netee of the teer. points v';) selected iron the first bend of the bese sc«ne (bl)t 
^^*b21*^b2p '■***’^****‘‘'^ *** coordinetes of the test points selected 

from the second bend of the bese scene (b2) • end t »y»< ) 1 represents the co 

ordinetes of the test points selected fr<» the aerged scene (a). These consti- 
tute ell the aeesured veriebler. Severel foxaules will be derived using the x 
coordinates. Similar forau!es exist for the y coordinates but will not be ex- 
plicitly derived. 


For the purposes of this aodel several ess\is|>tlons are aade. First it is 
assumed that, for e i, 

where is the true location of the test point and is the error of the 
htiaan selecting the points. Likewise » assume that 

*b2i " *i S2i 
Further* assuae that 


'ml 


*1 + + “t 


%fhere is the human error in selecting the point and o^ represents the mis- 
registration at that point. 


In addition, assuae that e^ 2 * ^ lndepend«it random 

2 2 2 

variables, normally distributed, with mean 0 and variances • o ^|^2 * °xm * 
and It is assumed that ^xa equal to o^. which will 

represent the hwsan error in selecting points. 
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Then consider the dlffsrsaess bstnera end psrtlculsr* 


S ^*bU~*b21^ 


1-1 n 

which can ba sioplifiad as follows: 
n 




* Sli ~*i “S2i^ 


i-1 n 

The assua^tion of independence yields 
n n 


E w . E 


®b21^ 


i-1 n i-1 n 

2 2 2 

" 5 V, + ® wo ■ 2o . 

xbl xb2 xh 

Thus ve get a siiqple formula for estimating the human error cogq>onenti 

o_J - 1/2 f; ~ %2i^^ 


xh 


i-1 


( 1 ) 


Now coiislder the difference between x^ and the average of x^|^ and x^2i* 


1-1 


( -,1 - 


4 I?! ( *1 ♦ * «t 


n 2 

c. + Oi - Sli - S2i ) 

2 2 


* " 

s 2 ( 

” 1-1 


\ 

I 
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Th« atmaptlon of lAdopondoiieo yloldit 


D 


E 




+ 1/*2 

i»i 



+ 1/4 



•V* »xh' 


2 

Ualng fonwl* (1) to roploco wo got: 



n 


n 


( 2 ) 


ThoM two foraulos ylold valuos for oiMi to toms of pixels, nio 
rns error of the control points Is repressed in Mters so sn expression of the 
Bisrsgis tret ion to asters is desirsble. 


"x ■ 


Slailsrly* a fonnals could be derived for a • Coablning o, end a * one gets 

y * y 


a value for the total ras error, o«: 


V"x^ * V 


A useful expression esn be deri^^ to deteraine the radius which contains 
90Z of the dats. Assiming a noraal distribution, the density function is 




2vo 


This esn be approxtoated as 


2 2 

2te' 


where a ■ 


0+0 

- » r . 
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IntagrAtiAg ov«r a elrela of radius X giuaa 


Q 




dx dy 



I - 

a r dr dO 

2to^ 



Solving thia for X glvaa 


X 



or 


R 




For Q *0.9 tha radiui bacoaaa 

O +0 

R 90 - (JL_X) 2.146 • 
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V. RB8ULTS AMD 00Mat»I(MI8 

TIm r«»itlt« of thio study srs listsd ia Tsbls 2. For tiM 12 sitss tsst^» 
ths 1LA8 slgoritla had an avaraga aargad scana arror (o^) of 31.6 aatara. In 
tha sight tost ottaa with coaparahla valuaa for tha BDim algorltla, tha avar> 
aga of ILAS 0 |^waxuas ««a 32.6, whlla tha avaraga of EDITCMl ^luaa waa 40.1. 
F>irthar, aach valua la lowar for tha ELAS data aato. Thus aa F-taat for 
olgnlficanca mis dona to trot whatlMr tha CAS waluaa wars a^lflcantly SMllar. 
Tha raaults ara llatad In Tabla 2. 

An aiwlysla waa dona to datamlna If any ralatltmahlp axlotad batwMO raa 
arrora and valuaa. ^wavar, this failad to ylald may at^stmtlal rasulta. 

Tha Spaarara's rank corralatlon coafflclant ba tw a an thaaa two varlablas waa 
0.159, a valua too low for olgnlficanca. 

Tha last analysis was to aaa if a significant dlffaranca axlatad for o 

X 

valuaa aaong tha four gmaral citagorlaa. Howavar, In all four catagorlaa 
thara was a algnlflcMt dlffarmca batwaan mom pair of valuaa. Thus dlffar- 
ancaa aaong catagorlaa ara laaa than diffarancaa within catagorlaa. 

In conclusion, una any stata that tha ERL scana-to-acma raglatratlon la 
significantly bat tar In flva of tha sight conparlamio. Thara ara two poaolbla 
raasona for tha dlffarMtcas. First, tha EC. aodal la a plMswlsa linaar aodal 
and tha EDITOR nodal is a cubic polynMlal aodal. Sacond, tha BKL progm ra- 
oaaplao using bUlnaar Intarpolatlon idilla tha OITOl sof twara ums a naaraat 
nalghbor raaaapllng. This sti^ did not Indlcata how auch ol' tha dlffarmea la 
attrlbut^la to aach factor. 
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TABLE 2. TEST RESULTS 


1 


ACCURACY AS 

SESSHENTS I 

Land Cover 
Test Areas 

ERL Control 
Point RMS 
Error (e^) 

Nuinber of 
Test Points 

NSTL/ERL SRS 

R90 

Model RMS 

(5m) 

Error 

F-’’est Results* 
(d.<f) ERL Vs. SRS 

NSTL/ERL 

SRS 

NSTL/ERL 

SRS 

Asrlculture 









Arizona 

34.9 

43 

43 

50.1 

65.1 

35.008 

43.521 

Significant 

South Dakota 

49.7 

52 

64 

28.9 

64.3 

26.916 

42.753 

Significant 

Missouri 

59.8 

58 


61.1 


40.401 



Ranse 









Kansas 

55.4 

64 

64 

55.9 

56.7 

37.603 

38.813 

Not Significant 

Idaho 

45.7 

63 

64 

48.2 

69.4 

31.810 

46.058 

Significant 

Oklahona 

49.4 

70 


42.0 


28.320 



Forest 









Colorado 

44.0 

63 

63 

54.4 

56.9 

36.875 

39.035 

Not Significant 

South Carolina 

49.4 

64 

64 

50.7 

49.4 

33.488 



Idaho 

49.6 

72 


39.0 


25.741 

32.575 

Significant 

Mixed 









Kentucky 

39.8 

62 

62 

49.3 

56.5 

32.894 

39.552 

Significant 

Horth Carolina 

39.0 

72 

64 

51.0 

58.3 

33.691 

38.447 

Not Significant 

Mississippi 

39.8 

71 


24.5 


16.236 



Averaiie 

46.4 



46.3 

59.6 

31.582 

40.094 



*®0* ^SRS * ^EBL • ® 


“a= »^S«S » •^WKL 


